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ABSTRACT 

To study the crucial range of Galactocentric distances between 12 and 16 kpc, where little 
information is available, we have obtained VI CCD imaging of Berkeley 20 and BVI CCD 
imaging of Berkeley 66 and Tombaugh 2, three distant, old open clusters. Using the synthetic 
colour magnitude diagram (CMD) technique with three types of evolutionary tracks of dif- 
ferent metallicities, we have determined age, distance, reddening and indicative metallicity 
of these systems. The CMD of Be 20 is best reproduced by stellar models with a metallicity 
about half of solar (Z=0.008 or 0.01), in perfect agreement with high resolution spectroscopic 
estimates. Its age r is between 5 and 6 Gyr from stellar models with overshooting and be- 
tween 4.3 and 4.5 Gyr from models without it. The distance modulus from the best fitting 
models is always (m-M)o=14.7 (corresponding to a Galactocentric radius of about 16 kpc), 
and the reddening E(B — V) ranges between 0.13 and 0.16. A slightly lower metallicity 
(Z~0.006) appears to be more appropriate for Be 66. This cluster is younger, t=3 Gyr, and 
closer, (m-M)o=13.3 (i.e., at 12 kpc from the Galactic centre), than Be 20, and suffers from 
high extinction, 1.2^ E(B — V) ^1.3, variable at the 2-3 per cent level. Finally, the results 
for To 2 indicate that it is an intermediate age cluster, with r about 1.4 Gyr or 1.6-1.8 Gyr 
for models without and with overshooting, respectively. The metallicity is about half of solar 
(Z=0.006 to 0.01), in agreement with spectroscopic determinations. The distance modulus is 
(m-M)o=14.5, implying a distance of about 14 kpc from the Galactic centre; the reddening 
E(B — V) is 0.31-0.4, depending on the model and metallicity, with a preferred value around 
0.34. 

Key words: Galaxy: disc - Hertzsprung-Russell (HR) diagram - open clusters and as- 
sociations: general - open clusters and associations: individual: Berkeley 20, Berkeley 66, 
Tomabugh 2 



1 INTRODUCTION 

Open clusters (OCs) are very good tracers of t he Galactic disc 
prope r ties, of its formation and evolution 
198 it iFriell 1 19951: 



(e.g. Panagia & Tosi 



Twarog, Ashman, & Antho nv-Twaroa 119971 : 



Freeman & B land-Hawtho rr l2002l) . In particular, OCs can be used 
to study the metallicity distribution in the disc and its possible evo- 
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lution with time. With the BOCCE (Bologna Op en Cluster Chem- 
ical E volution) project, described in detail by Bragaglia & Tosil 
fcOOfjl) . we are deriving precise and homogeneous ages, dis- 
tances, reddenings and chemical abundances for a large sam- 
ple of OCs. The final goal is to study the present status of the 
Galactic disc, its formation and evolution. Since the least known 
epochs of the disc evolution are the earliest ones, we put par- 
ticular attention to the study of old clusters. A dding the present 
three old systems t o those already describ ed in iBragaglia & Tosil 
'2006h, to Be 17 llBragaglia et alj [2006ah, Be 32 and King 11 



Tosi. Bragaglia & Cignonil l2007i), we have already examined 17 



clusters o l der th an 1 Gyr, out of the about 190 listed in the 
iDias et al.1 d2002h catalogue. 

As part of this project, we present here a photometric study 
of three clusters: (a) Berkeley 20 (I = 203. °48, b = -17.°37, 
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in the third Galactic quadrant); (b) Berkeley 66 (I = 139.° 43, 
b = 0.°22, in the second Galactic quadrant); and (c) Tombaugh 2 
(I = 232. °83, b = -6.°88, in the third Galactic quadrant). They 
are old, distant clusters, and their properties are important to un- 
derstand the nature of the outer Galactic disc. In fact, they all lie in 
the region between about 12 and 16 kpc from the Galactic centre, 
where the radial metallicity distribution seems to change its slope 
(see Sect.[6]l. Furthermore, To 2 has been connect ed with the Mono- 
ceros ring and the Canis Ma jor overdensity (e.g.. lFrinchabov et al .1 
l2004l ; lBellazzini et al]2004l see also Sect. l2.3t . 

Our paper is organized as follows: in Sect.[2]we give a short 
description of what is already available on these clusters in litera- 
ture, while in Sect. [3] we describe our data and the reduction pro- 
cedure. In Sect. [4] we present the CMDs and in Sect. [5] we detail 
the cluster parameter derivation. Finally, Sect. [6] is dedicated to a 
summary and discussion, in particular of the radial Galactocentric 
metallicity gradient and the importance of the three clusters in this 
context. 



low resolution spectroscopy were presented by IFriel et al. I J2002t) 
for nine stars; six of them appear to be cluster members and 
have (RV) = +70 ± 13 km s" 1 , [Fe/H]=-0.61 ± 0.14 dex. 
iFrinchabov et alj d2006l) observed 20 objects in Be 20 at interme- 
diate resolution and derived precise RVs; only five stars resulted 
members and their average RV is +75.7 ± 2.4 km s _1 . 

Abundan ces based on high resolution spectroscopy w ere first 
measured by lYong. Carney. & Teixera de Almeid3 d2005 ). They 
obse rved four stars, two in c ommon with IFriel et al. I d2002h . one 
with IFrinchabov et al.l d2006l) : their average RV is +78.9 ± 0.7 
km s - . Abundance analysis was possible only for two stars, 
located nea r the R GB tip, giving [Fe/H]~ -0.49 ± 0.06 dex. 
lYong et alj d2005h derived detailed abundances and discussed 
this cluster in the framework of G alactic abundance t rends and 
different origin/population. Finally, ISestito et alj d2008l) obtained 
FLAMES/UVES spectra of six stars in the field of Be 20 over the 
large WFI area; two of them are confirmed cluster members, and 
their analysis gives a metallicity [Fe/H]=— 0.3 (rms=0.02) dex. 



2 THE THREE CLUSTERS IN THE LITERATURE 

The three clusters have already been partly studied in the past and 
we briefly summarise here the available information. We have re- 
trieved (from the WEBDA3 or the original papers) the values of 
radial velocity (RV) for stars in our catalogues. They are presented 
in Tables [T][2] and [3] together with photometric data, coordinates, 
and identifications. The information on membership and metallic- 
ity (see next) for the clusters will be used in the present paper to 
help in the selection of the best-fitting synthetic CMDs. 



2.1 Be 20 

The first calib r ated p h 

iMacMinn et alj dl994l) . who obtained V, I data on a 5.1x5.1 
arcmin 2 field using the KPNO 2.1m telescope. Their CMD is 
well defined, but show s an apparent lack of red clump (RC) stars. 
IMacMinn etai] d 19941) derived, using isochrones, an age of 6 Gyr, 
[Fe/H]~ -0.23, (m - M) v ^ 15.0 and E(V - I) ~ 0.16. They 
deduced a radius of 1.5 arcmin and a mass of about 1000 Mq. They 
considered Be 20 worth of further interest also because of its large 
Galactocentric distance (Roc — 15.8 kpc) and unusual position 
belo w the Galactic plane (about 2 . 5 kpc) . 

iDurgapal. Pandev. & Mohanl J200ll) presented B, V, R, I data 
obtained at the 104-cm Naini Tal State Observatory over a 6x6 
arcmin 2 field. At variance with the other photometric works, they 
do think they see a horizontal branch in the innermost region of 
Be 20. Using two different sets of isochrones they concluded for 
an age of about 5 Gyr, E(B -V) = 0.10, Z=0.008 (i.e. [Fe/H]~ 
—0.3), (m — M)v — 15.1, a cluster radius of about 2.5 arcmin 
and aRpc — 17.1 kpc. 

iMomanv et al. liooH) showed results on B, V photometry ob- 
tained with the Wide Field Imager (WFI@2.2m ESO-Max Planck 
telescope) on a much larger field of view (about 30x30 arcmin 2 ). 
These observations were intended to produce catalogues over fields 
of view appropriate to the FLAMES fiber spectrograph and the au- 
thors did not really discuss the cluster properties. 

Radial velocities (see Table QJ and abundances based on 



2.2 Be 66 

For Be 66, the only available photometry reaching th e main se- 
quence Turn-Off (MSTO) is by IPhelps & Janesl dl996l) . They ob- 
served a 5.1x5.1 arcmin 2 field using the KPNO 2.1m telescope 
with the V, I filters, obtaining a well defined CMD. They derived 
the following parameters: age = 3.5 ±1.0 Gyr, —0.23 ^ [Fe/H] ^ 
0, E(V - I) = 1.60 ± 0.05, (m - M) v = 17.40 ± 0.20 (imply- 
ing a Galactocentric distance Rgc — 12.9 kpc), radius of 1.2-3.5 
arcmin and minumum mass of ~750 Mq; they also suggested the 
poss ibility of differential re ddening. 

IVillanovaetal.ld2005T) obtained high resolution spectra of two 
red clump stars with the HIRES spectrograph on the Keck I tele- 
scope. They have RVs —50.6 and —50.7 km s _1 (see Tablef2J, and 
seem to be cluster members. Abundance analysis was possible only 
for one star; Villanova et al. estimated [Fe/H]= —0.48 ± 0.24 and 
claimed that other elements have solar scaled ratios. 



2.3 To 2 

Although with significant discrepancies, all authors concur that 
To 2 is a d istant, old, rather m etal-poor open cluster. To 2 was dis- 
covered bv lTombaughl dl938l) and its first C MD, barely reaching th e 
main sequence turn-off, was published by lAdler & Janesl dl982l) . 
B,V,I images were collected with a variety of cameras at the 1 - 
m and 2.5-m telescopes in Las Campanas bv lKubiak etall d 19921) . 
who intended to find variable stars. Their CMDs are of very good 
quality (see also Sect. I3.2| l and they derive a distance of 6.3 ± 0.9 
kpc and an age of 4 Gyr, assuming E(V — I) = 0.4 and a metallic- 
ity one tenth of solar. lPhelps et al.ldl994h included To 2 in their list 
of old OCs, assigning it a value of SV = 1.5, i.e., an age of about 



2.5 Gyr (following the formula in I Janes & Phelpsl 



1994. No fur 



1 \http://www. univie. ac.at/ webda/webda. html see IMe rmilliod & PaunzenI 
120031) 



ther ph otometric catalogues are freely available, but lVillanova et all 
d2010h re-determined the parameters of To 2 on the basis of spec- 
troscopy (see below) and photometry, deriving an age of 2 Gyr, 
(m - M) v = 15.1, E(B - V) = 0.25, with [Fe/H]=-0.32 dex 
(see next). 

The metallicity of this cluster has b een d etermined using low- 
resolution spectroscopy by IFriel et alj d2002l) . who find [Fe/H] = 
—0.4 ± 0.09 using 12 stars; the corresponding RVs are given in 
Table [3] together with the ones by others, for a l l star s in com- 
mon with our photometric catalogue. IBrown et alj Jl996l) analysed 
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To 2, Be 20, and Be 66 3 



Table 1. Stars in our Be 20 photometric catalogue for which RVs have been published. FLAMES stands for the RVs from Sestito et al. (2007). ID in first 
column is our identificator; IDi, ID2 are in the WEBDA system; ID3 is taken directly from the Frinchaboy et al. (2006) paper. The precision of the velocities 
is 0.5, 10, 1, 2 km s" 1 for RV , RVi, RV 2 and RV 3 , respectively. 



ID 


V 


I 


RA(2000) 


DEC(2000) 


FLAMES 


Friel et a 


. 2002 


Yong 


et al. 2005 


Frinchaboy et al. 2006 


member 








hh:mm:ss 


dd:pp:ss 


RV 


IDi 


RVi 


ID 2 


RV 2 


11J3 


KV3 
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14.792 


13.325 


05:32:3 /.So 1 


■ c\f\ .11. fin £. 1 

+00:11:09.61 
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15.145 


13.749 


05:32:38.960 


+00: 11:20.33 


78.6 
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Y 
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14.565 


05:32:39.228 


1 f\f\ .1 A.11 f\ A 

+00:10:31.04 




10 
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10810 


63.7 


Y 
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15.752 


05:32:36.886 


1 f\f\ .1 1 . A r» AC\ 

+00: 1 1:49.49 




22 
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22 


78.9 


_ 


73.7 


Y 


1 85 


17.080 


15.928 


05:32:33.213 


+00:09:34. /9 




27 


1 5 






10617 


0.1 


N 
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16.841 


16.023 


05:32:42.328 


+00:12:29.39 




29 
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Y 
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17.446 


16.312 


05:32:37.014 


+00:12:00.54 




33 
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10741 


74.5 


Y 


281 


17.543 


16.442 


05:32:33.276 


+00:10:53.32 




37 


69 










Y 


780 


17.496 


16.492 


05:32:33.509 


+00:11:37.62 




39 


69 
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17.613 


16.534 


05:32:46.573 


+00:08:41.92 
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30.8 
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72 


17.040 
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+00:09:27.31 












10806 


35.3 
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05:32:46.420 
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85.4 














N? 


702 


16.187 


14.943 


05:32:36.774 


+00:11:04.84 


78.5 


1 1 












Y 


977 


15.447 


13.992 


05:32:50.067 


+00:16:16.59 


37.7 














N 



Table 2. RVs for the two stars (both members) in Be 66 analysed by Villanova et al. (2005). ID, magnitudes, and coordinates are from our work. 



ID 


IDwEBDA 


B 


V 


I 


RA(2000) 


DEC(20()0) 


RV 












hh:mm:ss 


dd:pp:ss 


km s — 1 


6521 


785 


20.375 


18.232 


15.596 


03:04:02.90 


+58:43:57.0 


-50.7 ±0.1 


5793 


934 


20.357 


18.217 


15.595 


03:04:06.41 


+58:43:31.0 


-50.6 ± 0.3 



high-resolution spectra of three stars, finding an average [Fe/H]= 
-0.4 ± 0.25 (with E(B-V) = 0.4), and rather normal elemental 
ratios (e.g., slightly enhanced a-elements, slightly deficient oxy- 
gen) although with a caveat on the uncertainties. 

This cluster could even be of extragalactic origin, and only 
recently accreted by the Milky Way: on the basi s of positional 
and kinematical arguments, iFrinchabov et alj J2004l) proposed that 
To 2 is part, together with other OCs, of the Mo noceros stream, 
also called GASS , Grea t Anticenter Stellar Stream l lNewberg et al .1 
120021 ; Ibata et al .1 120031) . poss ibly due to a dissolving, merging 
galaxy. iBellazzini et alj ( 120040 associated To 2 directly to the dis- 
solving dwarf galax y proposed as origin of the Monoceros stream 
dMartin et al.ll2004l) . i.e., to Canis Major (C Ma), whose reality ha s 
however been questioned by others (e.g.. iMomanv et all l2004h . 
To 2 has then attracted the attention of several groups and two pa- 
pers on its metallicity have recently appeared. 

IFrinchabov et alj J2008t) analysed high-resolution UVES and 
GIRAFFE VLT spectra of 40 stars, of which only about one half 
turned out to be members on the basis of the RVs. Their surpris- 
ing results is the presence of two sub-populations, both appar- 
ently member of the cluster, one more metal-poor, a-rich ([Fe/H= 
—0.28, [Ti/Fe]= +0.36, seven stars) and more centrally concen- 
trated, the other more metal-rich, with solar-scaled a elements 
([Fe/H= -0.06, [Ti/Fe]= +0.02, 11 stars) and more external. 
If confirmed, it would be the first case of chemical inhomogene- 
ity in OCs and maybe of multiple stellar generations. To 2 would 
in any case represent a very different case from globular clusters 



(GCs). There, inhomogeneities, initially regarded as "anomalies", 
have been found since a long time, but only in light elements, like 
C, N, O, Na, etc, whose abundance variations are (anti-)correlated 
dGratton. Sneden. & Ca rretta 200^). In contrast, the metallicity of 
GCs, described by [Fe/H ], is homogeneous to better than 10 per 
cent l lCarretta et al]|2009t) , with very few exceptions, like uj Cen- 
tauri, or M22. While the presence of multiple stellar generations 
in (probably all) GCs is presently widely accepted and i s based on 
both photometric and spectroscopic evidence (see, e.g.. lBragaglial 
l201fj|) . the mass of GCs is much higher than that of OCs, and 
mass is certainly of paramount importance in sh aping the destiny 
of a cl uster. To explain the peculiar case of To 2, IFrinchabov et alj 
d2008l) proposed several alternative solutions: two overlapping, or 
merged, clusters, multiple star formation periods, or the presence 
of both an open cluster and a stream, remnant of the same dwarf 
galaxy where also To 2 was born. 

This result has howeve r been questioned by part of the same 
team; IVillanova et alj d201(]|) presented results on GIRAFFE VLT 
spectra for 37 RGB and red clump stars (only 15 actually mem- 
bers, and only 13 with data good enough for the analysis). They 
found [Fe/H]= —0.31 ± 0.02 dex (rms=0.07 dex) and no sign of a 
bimodal distribution. All the other measured elements are also uni- 
form. They discuss the cause of this very discrepant result and the 
only reasonable explanation appears to be the better quality of the 
ne w data and the fac t that t heir spectra are in a redder region than 
the IFrinchabov et all d2008l) ones, hence more easily analysed for 
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Table 3. RVs for stars in To 2. ID, magnitudes, and coordinates are from our work. RVi: Brown et al 1996; RV2: Friel et al. 2002; RV3: Frinchaboy et al. 
2008; ID 4 and RV 4 : Villanova et al. (2010) (their Tables 1 and 5). 



ID 


IDwEBDA 


B 


V 


I 


RA (2000) 
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Table 3. (continuation) 
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relatively metal-rich and cool (i.e., with very crowded spectra) stars 
as the RGB ones in To 2. 

Whether CMa, and perhaps the Monoceros ring, are truly of 
extragalactic origin, or are instead due to disc warping and flaring, 
or spiral arms, or a combination of these phenomena, and whether 
To 2 is truly associated to these structures, or is "simply" a normal, 
Galactic-disc open cluster, the study of this part of the Galactic 
disc is important. Further observations, both photometric and spec- 
troscopic, are then welcome, in particular to settle the question of 
the possible chemical inhomogeneity, 



3 OUR DATA 

The three clusters were observed from two observatories and with 
three telescopes; a log of the observations is given in Table 4. The 
final photometric catalogues will be made available at the WEB DA. 
They will present the (B) Vic magnitudes with errors, pixel coor- 
dinates, and equatorial coordinates (tied to the GSC2 Catalogue us- 
ing CataPack, a software written by P. Montegriffo at the Bologna 
Observatory). 

3.1 Be 20 and Be 66 

Observations of the two clusters and the two control fields were 
obtained at the Telescopio Nazionale Galileo (TNG) on Canary 
Islands, during two nights in October and November 2000. We 
used DOLORES (Device Optimized for LOw RESolution), char- 
acterized by a field of view of 9.4x9.4 arcmin 2 and a scale of 
0.275 // /pixel. Since in these two runs we also observed B e 17 
jBragaglia et al.l2006 j) and NGC 6939 l lAndreuzzi et al.l2004l) . we 
refer to those publications for details on observation and reduc- 
tion procedures. Very brie fly, we used the package DAOPHOT-II 
dStetsonl 1987l : lDavisl 1994b to obtain PSF magnitudes, we corrected 
them to the same scale of aperture photometry and calibrat ed to the 
stand ard Johnson-Cousins system using Landolt's areas dLandoltl 
1 1992b . The calibration equations adopted for Be 20 and its control 
field observed on the same night, are: 

V = v - 0.0947 x (v + 1.1917 (rms = 0.014) 



Ic = i - 0.0060 x{v-i) + 0.7638 {rms = 0.025) 

The calibration equations adopted for Be 66 are: 
B = b + 0.0475 x (6 - v) + 1.4211 (rms = 0.012) 

V = v- 0.0959 x (v -i) + 1.2003 {rms = 0.012) 

Ic = i + 0.0422 x (v - i) + 0.7248 {rms = 0.012) 

The calibration equations adopted for the control field of 
Be 66 are: 

B — b + 0.0421 X (6 - v) + 1.4043 {rms = 0.012) 

V = v- 0.0947 x(v-i) + 1.1917 {rms = 0.014) 
Ic = i + 0.0017 x {v-i)+ 0.7589 {rms = 0.024), 

where b, v, and i, are the aperture corrected instrumental mag- 
nitudes after correction also for extinction and exposure time, and 
B, V, and Ic are the output magnitudes, calibrated to the Johnson- 
Cousins standard system. 

Finally, the completeness level of our photometry was derived 
with extensive artificial stars experiments, as in previous papers of 
this project. The completeness factors are shown in Table|5] From 
these experiments we also derived the errors typically associated to 
each magnitude level, that are used in Sect.|5]to build the synthetic 
CMDs. 

We compared our photometry for Be 20 with the ones of 
the two papers presenting the derivation of cluster parameters and 
readily availab le through the WEBD A, i.e., to the iMacMinn et all 
( 1 1994b and fhe lDurgapal et al.1 d200lb data. 

We counter-identified stars in the different catalogues and the 
results are presented in Fig. Q] In particular panels (a) and (b) in 
the figure show the comparison to the MacMinn et al. (1994) data, 
while panels (c) and (d) show the comparison to the Durgapal et 
al. (2001) data. For each star in common between the three cata- 
logs, V and I are the magnitudes in our photometry calibrated with 
the equations given in the te xt; Vaj and 1m are the corresponding 
magnitudes in the catalog o f lMacMinn etai] ( 1 1994b while Vd and 
Id are the magnitude s of the corresponding stars in the catalog of 
iDurgapal et aLU200lh . 
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Table 4. Log of observations for the clusters and the control fields 
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Figure 1. Comparison between our photometry and literature data for 
Be 20. Panels (a) and (b) show the comparison to the MacMinn et al. (1994) 
data; panels (c) and (d) show the comparison to the Durgapal et al. (2001) 
data. 
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For Be 66, the comparison was done to lPhelps & Janesl dl996t) . 
again only in the V and / filters, since they did not obtain B pho- 
tometry. The results are shown in Fig. [2] For each star in com- 
mon between the two catalogs, V and I are the magnitudes in 
our photometry calibrated with the equations given in the text. 
V p.j and Ipj are the corresponding magnitudes in the catalog of 
IPhelps & Janesl Jl996h - 

For both clusters, the differences between the different 
photometries have been computed retaining only well counter- 
identified stars (corresponding to the subsamples indicated in 
colour in the electronic version of the figures), with a difference 
in magnitude less than 0.15 mag, in absolute value. Differences are 
of the order of 0.02-0.03 mag in all cases. 

3.2 To 2 

Observations were obtained in La Silla, Chile, using two telescopes 
and three different instruments: a direct CCD camera (6.4x6.4 
arcmin 2 ) in 1995 and DFOSC (Danish Faint Object Spectrograph 
and Camera, 13x13 arcmin 2 ) in 2001, both mounted at the Dan- 
ish telescope, and SuSI2 (Super Seeing Imager 2, 5x5 arcmin 2 ), 
mounted at the New Tecnology Telescope (NTT) in 2002. We ob- 
served a field centred on the cluster and one about 50 arcmin away 
as a comparison to separate cluster and field stars. 

T he data reduction follow ed a standard procedure, see for in- 
stance |Pi Fabrizio et al] d2005h . Unfortunately, none of the nights 
turned out to be truly photometric, as we found out by com par- 
ison to the publis hed photometries by iKubiak et al .1 dl992h and 
IPhelps etal1dl994h . We decided to calibrate our photometry to the 
one by Kubiak et al., since it was the deepest one. No similar solu- 
tion has been possible for the comparison field; however, the pho- 
tometric accuracy obtained for it is enough for our goals (we only 
wish to separate cluster and field stars), also because we only have 
quite short exposures for this field. 

The results for the completeness tests are shown in Table [5] 
The final catalogue for To 2 contains 6073 stars. The limiting mag- 
nitude depends on the distance from the centre, since we obtained 
only short exposures with the widest field. In the following we will 
consider only the central 5x5 arcmin 2 region, observed with SuSI2 
on NTT, to derive the cluster properties. 



Figure 2. Comparison between our photometry and the one by Phelps & 

Janes (1996) for Be 66 * Be 20 - Fig.[3]shows the CMD obtained for the field centred 



4 THE COLOUR - MAGNITUDE DIAGRAMS 

• Be 20 - Fig. [3] shows the CMD obtained for the 
on Be 20 (left panel) and for the control field (right panel); the mid- 
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Table 5. Completeness level for the clusters and control fields of Be 20 and Be 66 and for the central field of To 2 as measured from the SuSI2 observations; 
mag is the calibrated magnitude (B, V, or Iq). 
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Figure 3. CMDs for Be 20 (left panel) and the control field (right panel). 
In the middle panel we indicate with different symbols member (blue open 
circles) and non member (orange crosses) stars, according to their RVs. The 
latter are listed in Table[T] together with the references. 



die panel indicates (with larger symbols) the stars for which infor- 
mation on membership is available thanks to the RVs. We clearly 
see the main sequence for about 4 magnitudes below the MSTO, 
located around V = 18.2, V — I = 0.7. We also identify the red 
clump with the few stars at V ~ 16.0, V — I ~ 1.2. Given the field 
stars distribution, stars above the MSTO may be attributed to a blue 
stragglers population or to fore/background contamination. Indeed, 
of the three stars with measured RV present in this region, two are 
field objects and one is a cluster member. More could be said only 
after decontamination, either statistical or (better) through actual 
measurement of the individual membership status via RV or proper 
motion. 
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Figure 4. Radial CMDs for Be 20 (upper panels) and the control field (lower 
panels). The stars within the 1', 2', 3' radii are 136, 269, 401 for the cluster 
and 19, 70, 153 for the control field, respectively. In the upper panels we 
show one of the best-fitting solutions (see Sect. l5.lt to guide the eye. 



Fig. [4] shows the radial distribution of stars in the cluster field, 
compared to an equal area in the control field; this is very useful 
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Figure 5. CMDs for Be 66 (left and middle panels) and the control field 
(right panel, shown only in the V, V — I plane). 
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Figure 6. CMDs obtained in different regions (500x500 pixels wide, i.e., 
2.2x2.2 arcmin 2 ) of Be 66. East is to the left and North to the top, and 
the cluster is located at the centre of the field of view. The isochrone (FST 
models with age 3 Gyr and metallicity Z=0.006) shown here represents one 
of the best solutions (see Sect. |5J and clearly indicates that differential red- 
dening is present, but only at a few per cent level. 



to better understand which are the true cluster sequences and the 
degree of field contamination (see also next Sect.). 

• Be 66 - The CMDs for Be 66 and the control field are pre- 
sented in Fig. [5] Notice the very red colours of the sequences, due 
to the high reddening, and the width of the evolutionary sequences. 
Ijanes & Phelpj {19941) suspected the presence of differential red- 
dening, a very reasonable possibility given its high value. We esti- 
mate that a variation of only a few per cent (see Sect l5.2b is suffi- 
cient to justify the observed spread. Fig. [6] shows a graphical rep- 
resentation of this: the CMDs obtained in sub-regions of the frame 
are all similar to each other and can be fit by the same isochrone 
(see Sect. |5.2| for the choice of the best one) once some difference 
in reddening is accepted. We notice in particular that the reddest 
CMDs correspond to the most eastern subregions in the two cen- 
tral rows of Fig. |6j while the least reddened CMDs are sparsely 
located: one in the bottom-left (south-east) panel and the other in 
the third top panel from left (north-east). This distribution suggests 
a clumpy/inhomogeneous extinction, rather than a systematic red- 
dening variation. 

In spite of the high dispersion and contamination, the main 
evolutionary phases, MS, MSTO, SGB, RGB, clump and possi- 
bly AGB are quite well defined, with the MSTO at V ~ 20.1, 
B — V ~ 1.6, and V — I ~ 2.1, and the clump - RGB intersection 
at V ~ 18.4, B - V ~ 2.2, and V - I ~ 2.7. 

• To 2 - The situation for To 2 is more complicated, since the 
cluster was observed with three different instruments with very dif- 
ferent field of view, and to very different depths (see Table 4). Fig. [7] 
shows the combination of all data for the central field in the upper 
row, left and middle panels, for the V,B — V and V, V — I CMDs, 
respectively. The entire control field is shown in the upper, right 
panel, only in V, V — I; we can immediately appreciate the differ- 
ent depths reached by the various instruments. 

To better separate the cluster CMD from the background, we plot 
in the lower panel of Fig. [7] only stars within a 2 arcmin distance 
from the cluster centre (left and middle panels) and within the same 
area in the control field (right panel). The evolutionary phases are 
much better recognizable, with the MSTO at V ~ 17.5, B — V — 
0.60, and V - 1 ~ 0.65, and the red clump at V ~ 16.2, B - V ~ 
1.12, and V-I ~ 1.33. 

Member and non member stars, according to their RVs (see Ta- 
bleO, are indicated with different symbols in the V, B — V CMD. 
This information is very useful to confirm the position of the RGB 
and red clump; unfortunately, RVs are not available for stars on the 
MS and the MSTO. However, as found bv lFrinchabov et all J2008L 
see their Sects. 3 and 6), even if the difference in RV between the 
cluster and the Galactic field stars is large enough to ensure a good 
decontamination, the same is not valid for the GASS/Monoceros 
component and some confusion may still be present. 

If we further restrict to the stars within 1 arcmin from the cluster 
centre, where the cluster stars dominate over the background, it is 
possible to see a clear indication of the presence of binaries, as 
shown in Fig.[8{a). The secondary MS is visible above and to the 
red of the single-stars MS; we use here B — I to have a larger 
baseline in colour. In Fig.[8tb) we indicate the two MSs; it is evident 
how the secondary, binary sequence complicates the definition of 
the position of the single-stars MSTO. 

The difference in magnitute between the MSTO and the red 
clump (8V) c an be used as an age indicator, after a suitable cal- 
ibrati on (e.g., IPhelps & Jane s|[l996l : iFrielll 19951 ; iBragaglia & Tosil 
l2006t) . This may be useful especially when dealing with large 
samples and/or a photometry that does not reach much below 
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Table 6. Values for SV, age in Gyr, Z, and [Fe/H] from high-resolution 
spectroscopy for the seven clusters analysed after Bragaglia & Tosi (2006). 
See text for the references. 



Cluster 


SV 


age 


Z 


[Fe/H] 


References 




Be 17 


2.9 


8.5 


0.008 


-0.10 


1,2 




Be 32 


2.6 


5.2 


0.008 


-0.29 


3.4 




Be 20 


2.2 


5.8 


0.008 


-0.30 


5,4 


King 11 


2.2 


4.5 


0.02 




3,- 




Be 66 


1.7 


3.8 


0.008 


-0.48 


5.6 


IE 


To 2 


1.3 


1.7 


0.008 


-0.31 


5,7 




NGC 3960 


0.8 


0.9 


0.02 


-0.12 


8, 8 





the MSTO . In Table [6] we indicate these values (measured as de- 
scribed in lBragaglia&Tbslll2006h for both the three OCs anal- 
ysed here and four others we published after iBragaglia & Tosil 
d2006h . We also give other relevant information, like age and 
metallicity (obtained with the synthetic CMD technique and the 
BBC tracks, see Sect. |5]( and [Fe/H]. The last column gives ref- 
erences for the parameters derived from photometry (Cols. 2-4, 
all by our group) and from high-resol ution spectroscopy (Co l. 



5) , respectively. The references a re: 1: IBragaglia et alJ 12UUbal) 
2: iFriel. Jacobson. & Pilachowskil ^20051) : 3: iTosi et alj jhoO/l) ; 4 



2006ah: 



Sestito et all j2008h; 5: this paper: 6: IVillanova et alj i2W3tl ; 7: 
Villanova et alJfeOlCh : 8: lBragaglia et alj J2006bl) . In Fig.l9lwe plot 
the values of SV versus age derived using stellar models for all OCs 
in the BOCCE sample. The relation is clearly not a simple one, 
as witnessed by the spread, only a pa rt of which is due to errors. 
Metall icity plays surely a role, see e.g. jTwarog & Anthonv-Twarod 
who decided to use a different indicator, combining differ- 
ences in magnitude and colour, to take it into account. We do not 
derive here a relation between 8V and age, deferring the task to 
when more OCs will be available on our scale. 
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5 CLUSTER PARAMETERS 

As fo r all the clusters of the BOCCE project (see Bra gaglia & Tosil 
120061 and references therein), we have derived age, distance and 
reddening following t he sy nthetic CMD method originally de- 
scribed by ITosi et all l ll99ll) . The best values of the parameters 
are found by selecting the cases providing synthetic CMDs with 
morphology, colours, number of stars in the various evolutionary 
phases and luminosity functions (LFs) in better agreement with the 
observational ones. To evaluate the goodness of the model predic- 
tions we quantitatively compare them with the observational LFs, 
stellar magnitude and colour distributions, and num ber of objects at 
the MS TO, the clump and th e RGB. As discussed bv lKalirai & Tosil 
d2004l) and lTosi et all J2007t) . where the luminosity and colour dis- 
tributions of each model were independently compared with the 
data using a Kolmogorov-Smirnov test and a \ 2 test, even sophis- 
ticated statistical procedures do not provide safer estimates of the 
cluster parameters. This is due to the background/foreground con- 
tamination and to the small number of objects usually measured 
in key evolutionary phases, such as the MSTO, the red clump and 
the SGB, fundamental in the identification of the cluster age. In 
the most favorable cases, the statistical tests confirm our choices of 
best synthetic models. In the cases of the three systems presented 
here, the parameter selection is even more uncertain, due to further 
decontamination problems related to the complicated and inhomo- 
geneous distribution of their fore/background interlopers. Yet, in 



Figure 7. V, B - V and V, V - / CMDs for To 2 and the comparison 
field. The upper panels show the entire field, while the lower panels show 
only stars within a 2 arcmin distance. We also indicate, only for the V, B — 
V diagram, member (blue open circles) and non member (orange crosses) 
stars, according to their RVs (see Table|3j. 
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Figure 8. (a) CMD in B, B — I of the central region of To 2 (within a radius 
of 1 arcmin) showing the binary sequence, (b) The same, but showing the 
MS ridge line and the same line shifted by 0.7 mag brighter. 
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Figure 9. SV versus age for the 23 OCs presently available in the BOCCE 
sample, 16 clusters from Bragaglia & Tosi (2006), seven from Table 6. The 
meaning of the different symbols is indicated in the figure. 



these unfavourable conditions the synthetic CMD approach is even 
safer (or less unsafe) than isochrone fitting than in standard cases, 
thanks to its capability to exploit all the available information on 
shape, population and position in the CMD of the various evolu- 
tionary phases. The method cannot provide strictly unique results, 
but allows to significantly reduce the range of acceptable parame- 
ters. 

In our procedure, the synthetic stars, extracted from the 
adopted stellar evolutionary tracks with a Monte Carlo approach 
taking into accout the adopted initial mass function (IMF) and star 
formation law, are attributed the photometric error derived from 
the artificial stars tests performed on the actual images. For all 
the BOCCE clusters we assume a Salpeter's IMF and a constant 
star formation lasting 5 Myr from the epoch of activation. The ex- 
tracted stars are retained in (or excluded from) the synthetic CMD 
according to the photometry completeness factors listed in Table[5] 
We have computed the synthetic CMDs both with and without tak- 
ing i nto account the possible contribution from unresolved binaries 
(see iBragaglia & Tosill2006l for a description of how binaries are 
included in the synthetic CMDs). Binary members are assumed to 
follow the same initial mass function (Salpeter's) as isolated stars 
and to have random mass ratio between primary and secondary 
components. In no case were the CMDs without binaries in agree- 
ment with the obervational ones. All the CMDs discussed here as- 
sume that 30% of the cluster measured stars are actually unresolved 
binaries. This fraction is consistent with what we find for the ma- 
jority of the BOCCE clusters . 

As usual, to test the effects of the adopted stellar evolution 
models on the derived parameters, we run the simulations with 
three different types of stellar tracks, with different assumptions 
for the treatment of convection, opacities and equation of state. The 
adopted models are listed in Table 7, where the corresponding ref- 
erences are also given, as well as the model metallicity and the 
information on their corresponding overshooting assumptions. 

To estimate the metallicity which better reproduces the pho- 
tometric properties of the cluster, we have created the synthetic 
CMDs adopting, for each type of stellar models, metallicities rang- 
ing from solar down to Z=0.004. We still assume as solar metal- 
licity models those with Z= 0.02, both for consiste ncy with the 
BOCCE previous studies (see lBragaglia & T osi 2006) and because 
they are the ones calibrated on the Sun by their authors. We con- 
sider only as indicative the metallicities obtained with our photo- 



Table 7. Stellar evolution models adopted for the synthetic CMDs. The FST 
models actually adopted here are an updated version of the published ones 
(Ventura, private communication). 



Set metallicity overshooting Reference 



BBC 


0.008 




yes 


Fagotto et al. 1994 


BBC 


0.004 




yes 


Fagotto et al. 1994 


BBC 


0.02 




yes 


Bressan et al. 1993 


FRA 


0.006 




no 


Dominguez et al. 1999 


FRA 


0.01 




no 


Dominguez et al. 1999 


FRA 


0.02 




no 


Dominguez et al. 1999 


FST 


0.006 


'/= 


=0.00,0.02,0,03 


Ventura et al. 1998 


FST 


0.01 


'/= 


=0.00,0.02,0,03 


Ventura et al. 1998 


FST 


0.02 


'/= 


=0.00,0.02,0,03 


Ventura et al. 1998 



metric studies and refer to high resolution spectroscopy for a safer 
determination of the chemical abundances. 

The transformations from the theoretical luminosity and effec- 
tive temperature to the Johnson-Cousins magnitudes and colours 
have been performed using the Bessell, Castelli & Pletz (1998 and 
private communication) conversion tables of the metallicity of the 
adopted models. We assume A v = (3.25 + 0.25 x (B - V) + 
0.05xE(B-V) )xE(B-V) and E(V - I) = 1.25 x E(B - V) x 
[1 + 0.06(5 - V)o + 0.0UE(B - V)] from Dean et al. (1978) 
for all sets of models. Using the same source for the conversion ta- 
bles for all models, we can be confident that the differences in the 
synthetic CMDs (and therefore in the cluster parameters) resulting 
from different stellar models must be fully ascribed to the intrinsic 
differences (input physics, opacities, etc.) of the models themselves 
and not to different photometric conversions. 

5.1 Be 20 

To minimize contamination without losing too many stars, we con- 
sidered as reference CMD the diagram of the stars located within 
2' from the cluster centre. As shown in Fig. [4] this is the region 
better compromising the two needs of removing contamination and 
having enough objects to define the key evolutionary phases. The 
CMD of this central circle is shown in the top panel of Fig. QJJJ It 
contains 269 stars, with 9 secure members (from the RV) nicely 
defining the red clump and the RGB, as shown in the central panel 
of Fig. [3] Since the control field of the same area contains 70 stars, 
we assumed that the cluster members within 2' are 200 and created 
the synthetic CMDs with this number of objects. 

Fig. [10] shows the synthetic CMDs in better agreement with 
the data for each set of tracks, selected by minimizing the dif- 
ferences between morphology, star counts, luminosity and colour 
functions in the observational and synthetic MS, MSTO, SGB, 
RGB and clump phases. The top panel of the Figure shows the ref- 
erence CMD, the panels in the central row show the best synthetic 
V, V - I CMDs obtained with the BBC, FRA, and FST models, 
overimposed to the CMD of the control field stars. The V, B — V 
synthetic CMDs corresponding to the same cases are shown in the 
bottom panels. 

As apparent from Fig.[4] field contamination affects mostly the 
CMD portion fainter than V ~18, but interlopers are present also 
at brighter magnitudes. Moreover, from that Figure and from the 
comparison of the top panel with the central row panels of Fig.fTOl 
the CMD of the control field stars turns out to be not quite the same 
as that of the cluster field contaminants. This circumstance, as well 
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Figure 10. Comparison between observational and synthetic CMDs for Be 20. Panel a shows the stars measured in V, and / in the central 2' radius region. 
Panels b, c, and d show the synthetic CMDs of best-fitting cases of each type of stellar model, overimposed to the CMD of the same area in the control field 
for a more direct comparison. Panels e, f, and g show the corresponding B — V CMDs. The displayed cases assume: Z=0.008, age=5.8 Gyr, (m-M)o=14.7, 
E(B - V)=0.13 for BBC, Z=0.01, age=4.3 Gyr, (m-M) =14.7, E(B - V)=0.16 for FRA, and Z=0.01, age=5.0 Gyr, (m-M) =14.7, E(B - V)=0.14 for FST. 



as the small number statistics, weakens the discrimination power 
of the LF comparison. In fact, all the reasonable cases present LFs 
similar to each other, such as those displayed in Fig.[TT] In this plot 
the dots show the LF of the 269 stars within 2' from the centre 
of Be 20, and the curves the LF of the (synthetic + field) stars of 
the CMDs in the central row panels of Fig. [10] The faint portions 
(dominated by field stars) of all the curves are in good agreement 
with the empirical LF, while the bright portions are consistent with 
it, but generally underestimated. We ascribe this underestimation 
of the number of relatively bright stars both to luminous interlop- 
ers and to the possible presence of BSS, not accounted for in the 
synthetic CMDs. 

Our photometry is only in two bands and, by itself, does not 
allow a reliable identification of the cluster metallicity. Models with 
solar metallicity or Z=0.004 do not reproduce the CMD morphol- 
ogy as well as models with intermediate metallicity, but cannot be 
excluded. Since the simultaneous agreement of the synthetic B — V 
and V — I with the corresponding observational col ours is a pow- 
erful metallicity indicator (see e.g. iTosi et alj|2007l for King 11), 
we have overcome our deficiency of information comparing the 
synthetic V, B — V diagrams with the corresponding CMD pub- 
lished bv lDurgapal et al.1 d200ll) . The self-consistency of this solu- 
tion is guaranteed by the agreement between our photometry and 
theirs, shown in Fig. [T] The V, B — V CMD of the stars from 



IPurgapaletail l bOOll) located within 2' from the cluster centre is 
shown in Fig. 1 121 

By comparing the synthetic diagrams based on all the stel- 
lar evolution sets and on various assumptions on the age, distance 
and reddening of Be 20 with our V , V — I CMD and with the 
V,B — V CMD from lDurgapal et al.1 J200ll) . we are able to distin- 
guish quite well all the cluster parameters. We find that only models 
with Z=0.008 or 0.01 can simultaneously reproduce the B — V and 
V — I observed colours. For higher metallicities B — V is always 
too red when V — I is fine, and, vice versa, for lower metallic- 
ities B — V is always too blue when V — I is fine. Hence, we 
assign to Be 20 a metallicity Z=0.009±0.001. This photometric es- 
timate of t he cluster metallicit y is in perfect agreement with that 
inferred bv lSestito et alj J2008I) from high-resolution spectroscopy 
([Fe/H]=-0.3). 

The age providing the best reproduction of the observed CMD 
morphology and stellar density in the different evolutionary phases 
depends on the assumptions of the stellar evolution tracks and is 
therefore slightly different from one set of models to the other. In 
spite of the field contamination and small number statistics, the 
clear identification from RVs of the SGB and clump stars allows 
us to infer the cluster age quite strictly. With the FRA models the 
age is in the small range 4.3 - 4.5 Gyr: older cases have the clump 
too bright, and younger cases have it too faint. With overshooting 
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Figure 11. Comparison between observational (dots) and synthetic (curves) 
LFs of the stars within 2' from the centre of Be 20. The curves include the 
200 synthetic stars and the 70 stars of the control field, and correspond to 
the cases shown in Fig. 1101 solid for the BBC models, dotted for the FRA 
and dashed for the FST ones. 
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Figure 12. V, B — V CMD of stars within 2' from the centre of Be 20, 
from Durgapal et al (2001). 



models, the age is obviously older, but only by less than 25% be- 
cause in these low mass stars overshooting is not very effective. 
With the FST models the best-fitting ages are between 5 and 5.5 
Gyr, while with the BBC models they are between 5.5 and 5.8 Gyr. 

Reddening and distance modulus are also well identified 
by best-fitting magnitudes and colours of the main evolutionary 
phases. Particularly striking is the circumstance that all the mod- 
els in agreement with the data require (m-M)o=14.7. E(B — V) is 
more sensitive to the metallicity and the details of the stellar mod- 
els, and we find it to range between 0.13 (for all the BBC-0.008 
good cases) and 0.16 (corresponding to the FRA-0.01 case with 
age 4.3 Gyr). 

In summary, thanks to the cluster membership provided by the 
radial velocities, age, distance, reddenning and metallicity of Be 20 
are very well identified. 

5.2 Be 66 

The situation for Be 66 is apparently more complicated, both be- 
cause of its high extinction and of the lack of information on mem- 
bership. The fraction of fore/background contaminating objects 
doesn't significantly vary from the cluster centre to the periphery 
and we have therefore chosen to simulate the whole CMD. Since 



the stars measured in all the three B, V and I bands are 2362 in our 
cluster field and 1023 in the control field, the synthetic CMDs have 
been created with 1339 objects, assuming the photometric errors 
and the completeness factors described in Sect. [3] 

The major problem encountered by all models is the width 
of the evolutionary sequences, both in colour and in magnitude. 
In no way are we able to create synthetic diagrams with the ob- 
served width by including only photometric errors and binary sys- 
tems. However, the synthetic sequences become properly thick if 
we assume variable amounts of reddening. We find that a SE(B — 
V) =±0.02 is sufficient to account for the observed spread. Given 
the high reddening affecting Be 66, this variation corresponds to 
only 2-3 per cent and looks rather likely. 

The widths of the MSTO, SGB and clump could also make it 
more difficult to precisely define the cluster parameters. However, 
it turns out that the independent constraints from the colour, mag- 
nitude, morphology, stellar density of the different evolutionary se- 
quences and luminosity function do allow to significantly reduce 
the range of possible values. The distance modulus, for instance, 
turns out to be strikingly stable in all the models in acceptable 
agreement with the data: in spite of its dependence on age, red- 
dening and metallicity, we find 13.2 ^ (m — M)o^l3.5 in all the 
acceptable cases, and (m — M)o=13.3 in all the best cases. 

To identify the cluster metallicity, the clue is the simultane- 
ous consistency of both the predicted B — V and V — I colours 
with the observed ones. All the models with solar metallicity and 
right B — V predict too blue V — I, so we can therefore ex- 
clude Z=0.02. For lower metallicity, the different sets of stellar 
tracks have the following responses: at Z=0.01, the FST models 
still show the B — V inconsistency with V — I of the solar 
models, although at a lower level, while the FRA models cannot 
be excluded. The BBC models with Z=0.008 have self-consistent 
colours and acceptable CMD and LF properties. The FST mod- 
els with Z=0.006 perfectly match the observational CMD and LF, 
and the FRA models with Z=0.006 also have self-consistent colours 
and accetable CMDs and LFs. Also the BBC models with Z=0.004 
give self-consistent colours and acceptable CMDs and LFs. Within 
the uncertainties, all these results indicate that the metallicity of 
Be 66 is most likely Z=0.006±0.002. This value corresponds to 
[Fe/H]~ —0.45, lower than the range proposed bv lPhelps & Janes! 
i 19961) but in perfect agreement wi th the clump star abundance 
measured by IVillanova et al.l d2005h from high resolution spec- 
troscopy. 

The reddening resulting from the synthetic CMDs depends 
both on the adopted models and metallicities. In the cases in better 
agreement with the data, we find it to range between E(B — V) = 
1.22 ± 0.02 and 1.30 ± 0.02, where 0.02 is the variation nec- 
essary to reproduce the observed spread of the evolutionary se- 
quences. In the best-fitting case (FST with Z=0.006 and r=3 Gyr) 
E{B — V) =1.26 ±0.02. 

The age depends on the various assumptions as well. How- 
ever, in spite of the width of the evolutionary sequences and the 
uncertainties on reddening and metallicity, we find it to be rather 
well determined for each kind of stellar models, thanks to the well 
defined shape and stellar density of key phases, such as MSTO, 
SGB, RGB and clump. With the FRA models without overshoot- 
ing we find 2.7 < r/Gyr <3.0, both for Z=0.006 and Z=0.01, 
with the preferred value between 2.7 and 2.8 Gyr. The models with 
overshooting provide 3.5 ^ r/Gyr ^ 4.0 with the BBC models 
with either Z=0.004 and Z=0.008 (preferred value 3.8 Gyr), and 
3.0 ^ r/Gyr ^ 3.5 with the FST models with intermediate over- 
shooting (preferred value 3.0 Gyr). 
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Figure 13. Comparison between observational and synthetic CMDs for Be 66. Panels a and b show the stars measured in B, V, and I in our cluster field. 
Panels c, d and e show the synthetic V, B — V CMDs of best-fitting cases of each type of stellar model, overimposed to the CMD of the same area in the 
control field for a more direct comparison. Panels f, g and h show the corresponding V, V — I CMDs. The displayed cases assume: Z=0.008, age=3.8 Gyr, 
(m-M) =13.3, E(B - V)=l.20, 1.22, 1.24 for BBC; Z=0.006, age=2.7 Gyr, (m-M)o=13.5, E(B - V)=1.28, 1.30, 1.32 for FRA, and Z=0.006, age=3.0 Gyr, 
(m-M) =13.3, E(B - V)=1.24, 1.26, 1.28 for FST. 




Figure 14. Comparison between observational (dots) and synthetic (curves) 
LFs of the stars in Be 66. The curves include the 1339 synthetic stars and 
the 1023 stars of the control field, and correspond to the cases shown in 
Fig. 1131 solid for the BBC models, dotted for the FRA and dashed for the 
FST ones. 



Fig[T3]shows the synthetic CMDs in better agreement with the 
V, B — V and V, V — I diagrams for each set of tracks, selected 
by minimizing the differences between morphology, star counts, 
luminosity and colour functions in the observational and synthetic 



MS, MSTO, SGB, RGB and clump phases. The two top panels of 
the Figure show the observational CMDs, the bottom panels show 
the best synthetic V,V—I CMDs obtained with the BBC, FRA and 
FST models, while the panels in the central row show the V,B — V 
CMDs of the same cases as in the bottom row, overimposed to the 
CMD of the control field stars. 

Fig|14l compares the cluster LF with the LFs from the syn- 
thetic cases shown in Fig.[T3] The three models predict LFs so sim- 
ilar to each other that the three curves overlap almost completely. 
They all agree very well with the data, except at the bright end, 
where field contamination dominates and the difference between 
the foreground stars in the cluster field and in the control field is 
apparent. 

For Be 66, the FST models with Z= 0.006, age=3 Gyr, E(B - 
V) =1.26±0.02 and (m — M)o=13.3 are by far the ones in better 
agreement with the data. 

5.3 To 2 

As discussed in the previous Sections, To 2 presents a fairly com- 
plicated CMDs, strongly affected by (probably multiple) contam- 
ination. We did observe a nearby region as control field to eval- 
uate the contamination, but the comparison of its CMD with that 
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of To 2 clearly shows that the two populations are rather different 
from each other (see Fig[7j. Moreover, we didn't have the oppor- 
tunity to acquire images of the external field in all the filters, and 
B is missing. Hence, the control field is of little help to decontam- 
inate To 2. From the morphology of To 2's CMDs with increas- 
ing distance from the cluster centre it is apparent that in the circle 
within 2 arcmin the contamination is still rather significant (at least 
15%, according to the number of stars in an equal area region of 
the control field). To minimize it, we thus decided to adopt as ref- 
erence cluster CMD that of the central region of 1 arcmin radius 
(see Fig[8}. It contains 1 147 objects with measured V and I magni- 
tudes, and 846 with B, V and I. Since control field portions of equal 
size contain 100 stars with measured V and I, we can assume a 9 
percent (lower limit to) contamination. In other words, of the 846 
objects with measured BVI, 770 can be considered To 2 members. 
The synthetic BVI CMDs have been therefore computed assuming 
this number of objects. The B-V and V-I CMDs of the 846 BVI 
objects in the central region of To 2 are displayed in the top panels 
ofFigE] 

Also in this case, we have tried to blindly identify the metallic- 
ity, without biassing our search on the basis of literature estimates, 
for homogeneity with the BOCCE project procedure. All the mod- 
els clearly exclude a solar metallicity: whatever the adopted tracks, 
Z=0.02 always lead to CMDs inconsistent with the data both be- 
cause of the MS shape and of the colours. If we assume a redden- 
ing allowing to reproduce the B-V colours, then V-I is systemati- 
cally too blue, and vice versa. Moreover, the acceptable reddenings 
are systematically lower than in literature. While rejecting a so- 
lar metallicity is straightforward, finding the best fitting one is ex- 
tremely difficult, because contamination affects key features before 
and after the MSTO. Since a good identification of the metallicity 
is necessary to identify the reddening, our analysis leaves the latter 
rather uncertain too. Vice versa, the distance modulus of To 2 turns 
out only moderately sensible to the metallicity choice and between 
14.4 and 14.7 for any viable model. 

With the BBC models, the available metallicity most appro- 
priate for the CMD of To 2 is Z=0.008, whose models allow to 
self-consistently reproduce the observed B-V and V-I colours of 
all the evolutionary phases. However, the corresponding upper MS 
and MSTO regions never have exactly the same morphology as 
the observed ones. The best CMD with these tracks assumes age 
t = 1.6 Gyr, reddening E(B-V) = 0.34 and distance modulus (m- 
M)o=14.5. Its CMDs are plotted in panels c and d of Fig|T5] Panel 
f shows only the 770 synthetic stars, while panel c contains 846 
objects like the empirical CMDs of the top panels, since it includes 
also the 75 objects falling in an equal area portion of the control 
field. The comparison of panels b and c emphasizes the difference 
between the control field and the cluster contaminating stars. BBC 
models with Z=0.004 also lead to an acceptable self-consistency 
between B-V and V-I colours, but systematically have RGBs red- 
der and MSs straighter than observed. For this reason we consider 
them less appropriate for To 2 than the Z=0.008 ones. 

For the other types of models, both the FRA and the FST 
metallicities below solar are Z=0.01 and Z=0.006, and for both it 
is impossible to significantly discriminate between them. With the 
FRA models, both metallicities allow for self-consistent colours, 
but both provide RGB always quite redder than observed. To shrink 
the subgiant branch and keep the RGB within the observed colour 
range, one should increase the age, but this inevitably implies 
clumps brighter than observed. Generally speaking none of the 
FRA models reproduces adequately the properties of To 2's CMD. 
In panels d and g of Fig[l5] we simply show the case in smaller 



disagreement, with no claim of actual consistency. It corresponds 
to Z=0.006, age 1.4 Gyr, E(B-V)=0.40 and distance modulus (m- 
M) =14.5. 

With the FST models, discerning the right metallicity is 
equally difficult, because Z=0.006 leads to the best self-consistency 
between B-V and V-I, but predicts RGBs slightly redder than ob- 
served, while Z=0.01 has V-I slightly too blue when B-V is correct 
(a typical signature of a metallicity overestimate), but appropriate 
RGBs. At variance with the FRA models, the FST ones provide 
however a good agreement with the empirical CMD, actually better 
than with any of the other models. Based on the overall properties 
of the CMD (curvature of the MS, shape of the MSTO and subgiant 
branch, colour and slope of the RGB, morphology of the clump, 
number of stars in each of these phases) we ended up preferring the 
Z=0.01 models. Most likely, the right metallicity (at least within 
the framework of the FST models) is slightly below this value. In 
panels e and f we have chosen to show one of the best cases among 
the FST ones: age 1.7 Gyr, E(B-V)=0.31 and distance modulus (m- 
M)o=14.7. Here the choice is rather subjective: an age of 1.6 Gyr 
(with same metallicity and distance modulus and E(B-V)=0.32), or 
an age of 1.8 Gyr (with same metallicity, E(B-V)=0.30 and distance 
modulus (m-M)o=14.5) would have been good as well. 

Since the FST models come in three overshooting flavours 
(see Table 7), we have tested them all. As often found for the 
BOCCE clusters, the FST models with overshooting reproduce the 
CMD of To 2 much better than those without. In this case, we found 
no significant difference in the quality of the fit from models with 
the highest or intermediate overshooting. 

We have tested various fractions of unresolved binaries and 
found that while their presence is clearly required to reproduce the 
secondary MS seen (see Fig[8]l at the right of the main MS both in 
B-V and in V-I (hence a binary fraction of must be rejected), we 
don't need to invoke particularly high percentages. A 60% fraction 
would definitely lead to an excessive blend of the binary and single 
star MSs and is thus beyond the acceptable value. We conclude that 
a fraction around 30% provides a good representation of the data. 

The difficulty in properly characterizing the field contamina- 
tion has the further consequence of reducing the selection power of 
the comparison between synthetic and empirical LF. In Fig|16lthe 
LFs of the three models plotted in panels c, d and e of Fig|15|are 
plotted as lines, while the empirical LF of the 846 stars of the cen- 
tral l'of To 2 is represented by dots with error bars. The synthetic 
LFs fit reasonably well the portions where the cluster dominates, 
but at both the faint and the bright ends, where contamination dom- 
inates and the control field is not representative of that population, 
all the LFs fail, inevitably unable to predict the number of interlop- 
ers. Notice however that in the central parts of the plot, the agree- 
ment is fairly good. In particular the bump and dip corresponding 
to the clump are well fitted, showing that both age and distance 
modulus are appropriate. 

To summarize, we conclude that within the high uncertainties, 
To 2 turns out to be a cluster of intermediate age (1.4 Gyr for mod- 
els without overshooting, 1.6-1.8 Gyr for models with overshoot- 
ing), with distance modulus 14.5-14.7, and with reddening between 
0.30 and 0.40 depending on the chosen metallicity, but most likely 
within 0.31-34 (i.e. for 0.008«SZsS0.01). 

It is interesting to notice that, in spite of the uncertainty in 
the identification of the best met allicity, we don't find an y need of 
a metallicity spread as found by iFrinchabov et all d2008l) . Assum- 
ing a metallicity or reddening dispersion in our synthetic CMDs 
would introduce a spread in the evolutionary sequences larger than 
observed. Notice that only their metal-poor population falls in the 
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Figure 15. Comparison between observational and synthetic CMDs for To 2. Panels a and b show the stars measured in B, V, and / in the central 1' radius 
region. Panels c, d, and e show the synthetic CMDs of best-fitting cases of each type of stellar model, overimposed to the CMD of the same area in the 
control field for a more direct comparison. Panels f, g, and h show only the synthetic stars of the corresponding B — V CMDs. The displayed cases assume: 
Z=0.008, age=1.6 Gyr, (m-M) =14.5, E(B - 10=0.34 for BBC, 2=0.006, age=1.4 Gyr, (m-M)o=14.5, E(B - V)=0.0.40 for FRA, andZ=0.01, age=1.7 Gyr, 
(m-M) =14.7, E(B - V)=0.31 for FST. 



central 1' -radius circle, so our result of a single population of stars 
with a single metallicity 0.008^Z^0.01 is in perfect agreement 
with theirs, as are also reddening and distance modulus. Since in 
the CMD of the region with 2'-radius we already see significant 
contamination, and since the most central metal-rich star is at least 
1.5'from To2's centre (see their Fig. 2), of the various alternatives 
examined by Frinchaboy et al. we favour the interpretation of the 
metal richer populations in terms of an interloper. The more recent 
findin g of absence of any metallicity dispersion dVillanova et al.l 
l20ld) is even more easily in agreement with our study. The metal- 
licity is similar, as are the distance, reddening, and age. 
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6 SUMMARY AND DISCUSSION 

We compared the CMDs of the three old, distant OCs Be 20, Be 66, 
and To 2 to synthetic ones based on three different sets of evolution- 
ary tracks and determined the clusters parameters. Table [8] gives a 
summary of the derived parameters for the three OCs. We use here 
the values based on the BBC tracks, even when they do not offer 
the best-fitting solution (e.g., for Be 66), in analogy to what we did 



Figure 16. Comparison between observational (dots) and synthetic (curves) 
LFs of the stars within 1' from the centre of To 2. The curves include the 
770 synthetic stars and the 76 stars of the control field, and correspond to 
the cases shown in Fig. 1151 solid for the BBC models, dotted for the FRA 
and dashed for the FST ones. 

in lBragaglia & Tosll d2006h to obtain a homogeneous ranking on a 
single scale. 

• All three clusters are less metal-rich than the Sun, with 
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Table 8. Summary of the parameters for the three clusters. We use the BBC 
results, for homogeneity in ranking and immediate compatibility with the 16 
OCs presented in Table 1 of Bragaglia & Tosi (2006); see text for the values 
derived using the other tracks. Note that Be 66 has differential reddening. 



oc 


age 


(m - M) 


d© 


Rac 


E(B - V) 


Z 




(Gyr) 




(kpc) 


(kpc) 






Be 20 


5.8 


14.7 


8.71 


16.0 


0.13 


0.008 


Be 66 


3.8 


13.3 


4.57 


12.0 


1.22 


0.008 


To 2 


1.6 


14.5 


7.95 


14.3 


0.34 


0.008 



best=fitting solutions of Z=0.008-0.01 (depending on the set of 
tracks) for Be 20 and To 2, and Z=0.006 for Be 66. These rather 
precise values could be obtained because we considered the simul- 
taneous good fit of both the V - I and B — V CMDs. 

• The reddenin gs slightly depend on th e tracks metallicity. For 
Be 66, we confirm Iphelps & Janesl dl996h 's suggestion of a prob- 
able differential reddening, of the order of 2-3 per cent, a rather 
likely occurrence given the high extinction level. 

• For all clusters a binary fraction of about 30 per cent seems 
necessary to well reproduce the width of the observed sequences. 

• The derived ages depend on the treatment of convection 
adopted in the evolutionary tracks, with the usual lower values 
found for tracks without overshooting. 

• Our results on age, distance and reddening do not significantly 
differ from the most recently published ones for Be 20 and To 2. We 
find Be 66 sligh t ly you nger, closer, and definitely metal poorer than 
Phelps & Janes ( 1996), w hile w e agree with them on the reddening 
and with lVillanova et al.l d2005l) on the rather low metallicity. 

• As already anticipated in Sect. [I] these OCs have Galactocen- 
tric radii of about 12 to 16 kpc and are then useful to constrain 
the properties of the outer disc, in particular in what seems to be a 
transition region for the metallicity distribution (see below). 

The Galactocentric metallicity gradient - Even if the number of 
known old OCs h as steadily grown in the last years, they still are a 
minority, since the lDias et al. d2002l) catalogue contains about 2000 
objects, and only about 190 are older than 1 Gyr (see Fig. 1 17b - 
An even smaller number of clusters has the metallicity determined 
using high-resolution spectroscopy (see below). In the BOCCE 
project we have especially targeted old clusters, hence our OCs rep- 
resent a fair sample of the old cluster population, about 10 per cent. 
In Fig. [T7] we also show for comparison the same histogram for 
our sampleQ and indicate the three clusters analysed here in the 
enlargement. 

The metallicity distribution in the disc is an important ingre- 
dient of chemical evolutionary models, and offers essential infor- 
mation on the formation and evolution of the disc, especially when 
also its possible evolution with time is considered. We may ob- 
tain the present-day metallicity distribution from 0,B stars and H 



2 This comparison is only indicative, since age determinations are not 
homogeneous between our work and the catalogue (and inside the cata- 
logue itself). There are differences, some times large , betw een ages mea- 
sured by different authors. For instance, in lDias et alj |2002j), a) one of the 
two clusters in the last bin -age 10 to 11 Gyr- is Be 17 that, according to 
iBragaglia et"aT1 {2006a ) is instead less than 9 Gyr old, b) King 1 1 is gi ven 
an age of 1.1 Gyr, while we derived an age > 4 Gyr (Tosi et al. 2 0071) . in 
line with other determinations, and c) NGC 6791 is attributed an age, 4.4 
Gy r, younger than fo und in all recent papers, where its age is 7-9 Gyr (see, 
e.g, lKingetalj2005l) . 
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Figure 17. Age distribution for the 1 86 OCs older than 1 Gyr, according 
to the 2010 update of the Dias et al. (2002) catalogue. The filled histogram 
represents the 17 clusters in the BOCCE sample within the same age limit 
(14 already published, three in the present paper). The inset shows an en- 
largement to better appreciate the -still small- number of very old objects, 
and the three new OCs are indicated by arrows. 



II regions (see, e.g.. iRudolph et al ] |2006h . Another possibility ire 



Cepheids (e.g.. lAndrievskv et alj2004f) . for which the distance and 
age can be determined with high precision. Planetary nebulae (PNe) 
are another tracer, and they are in princ iple able to reach further 
in the past (see, e.g., the recent paper bv lStanghellini & Havwoodl 
120101) ; however, there still are difficulties in assigning individual 
ages and distances. This is in general true for all field stars; the 
most accurate parallaxes to date, those measured by the Hipparcos 
space mission, are available only for the Sun's vicinity. Ground- 
based catalogues can reach farther away, but at the price of accu- 
racy. HST parallaxes can be very precise, but only a handful of 
targets has been observed. To obtain precise distances for hundred 
millions stars in the whole Galaxy we will have to wait about 10 
years, for the completion of the Gai^f] satellite survey. 

Furthermore, all stars are subject to orbit migration (see, e.g., 
Sellwood & Binn ev 2002; Roskar et al.ll2008l ; ISch6nrich & Binnevl 



2009L in the recent years); they move away from their birthplace 



even several kpc, thus complicating the study of the metallicity dis- 
tribution: at any given Rgc we may find stars born there, or in an 
inner (i.e., in general more metal-rich) or in an outer (i.e., in general 
more metal-poor) region of the disc. 

OCs are less subject to all these problems: their distances and 
ages can be measured with sufficient precision using the stellar 
mode ls (e.g., as we do in our BOCCE project, see lBragaglia & TosH 
120061) . their abund ances can be determin ed from several or many 
member stars (e.g.. IBragaglia et ak I l2008l) . and they do n ot appear 
to suff er from orbit migration. For the last point, see e.g JWu et al .1 
( l2009h . who computed the or bits of about 500 OCs, using infor- 
mation from lDias et alj d2002l) . When they derive the slope of the 
metallicity gradient using Rgc or the apogalacticon radius, they 
do not find any significant difference. Apart from a few exceptions, 



See http://www.rssd. esa.int/index.php?project=GAIA fcpage=index. 
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it appears safe to use OCs and their present-day positions, to define 
the metallicity distribution now and in the past. This is part of our 
goals, and we are building the BOCCE sample choosing OCs which 
cover the whole distribution of clusters' properties. Of course, the 
situation will further improve when it will be possible to compute 
orbits for the whole family of OCs, using the improved distances, 
proper motions, and radial velocities produced by Gaia. 

We show in Fig. 1181 a) the radial metall icity gradi e nt de- 
scribed by OCs using the information in the iDias et al] d2002h 
catalogue. We computed Roc using 8 kpc as the Sun dis- 
tance from the Galactic centre, and the tabulated distances 
from the Sun. Values for [Fe/H] come from photometry (CMD, 
narrow-band, indices), low resolution, and high resolution spec- 
troscopy; individual references can be found at the Dias web- 
page (http://www.astro.iag.usp.br/~wilton/!. We display with filled 
(red) symbols the postion of Be 20, Be 66, and To 2, using the 
tabulated values. The metallicity tends to decline from the cen- 
tre to the outskirts of the disc. If we fit the distribution with 
a single relation, we find a decline rate of —0.04 dex kpc -1 , 
as indica ted in the figure. However, recent studies of open clus- 
ters (e.gjYong et alJl2005l:ICarraro et alj|2007l ; ISestito et alj|2008l ; 
iFriel. Jacobson. & Pilachowskil 20 id) suggest that a better repre- 
sentation is a gradient in the inner region (^12 kpc from the cen- 
tre) followed by an almost flat value therefore. We then divided the 
sample of OCs accordingly, computing a (steeper) inner slope of 
—0.07 dex kpc - . The OCs in the external part have in this case 
an average metallicity of =-0.35 dex. While the number of OCs 
in this sample is large (177 objects), they tend to lie mostly in the 
Sun's vicinity, with only about 15 per cent of them with Roc > 12 
kpc, and only a handful in the outermost disc regions. Furthermore, 
the sample is completely inhomogeneous. 

In principle, choosing only [Fe/H] measured using high reso- 
lution spectroscopy produces more solid results. We have searched 
the literature and retrieved the metallicity, distance, and age of 
about 70 clusters. A good fraction of them, about 30 per cent, are 
also in the BOCCE sample and are indicated with filled symbols 
in the figure; we show the three clusters of the present paper with 
larger symbols, this time using the parameters derived in the present 
paper and the [Fe/H] values indicated in Table [6] We restricted our 
search to clusters older than 0. 1 Gyr and selected in each case ei- 
ther the most recent determination or our measures, when present. 
The resulting metallicity distribution is shown in Fig. 118( b). The 
slopes derived are (almost) the same as in the previous case. In this 
case, however, there is a better balance between inner and outer 
clusters, and the individual values are more precise (but again, with 
a caveat on the inhomogeneity of sources, which produces system- 
atics). Be 66, To 2, and Be 20, in order of increasing distance from 
the centre, are an important addition, since they are about one fifth 
of the entire outer-disc sample. 

This second sample of OCs can also be used to study the be- 
haviour of the gradient with time. In Fig.[T9]we divide the sample 
in three parts and compute the slope of the gradients (only in the in- 
ner 12 kpc) again, with the three clusters hightlighted. While young 
and intermediate age clusters seem to share the same behaviour, 
the slope of the gradient seems to have been steeper in the past (as 
found also for PNe). This has to be accounted for by any chemical 
evolution model. Note however the paucity of very old clusters at 
the transition between a decreasing and flat distribution of metallic- 
ities; this calls for new additions to the well-studied cluster sample. 

We plan to repeat these exercises using only clusters in the 
BOCCE sample, with ages, distances, and [Fe/H] all derived on the 
same scale. The importance of Be 20, Be 66, and To 2 lies both 
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Figure 18. The metallcity gradient as defined by OCs. (a) Rqc an d [Fe/H] 
are taken from Dias et al. (2002). The two lines represent the fit to the 
data using a single slope (in blue) or two, within 12 kpc from the Galactic 
centre and farther than this limit (in red), see text. The three clusters of the 
present paper are indicated by filled (red) squares. (b)The same, but using 
only 72 OCs older than 0.1 Gyr and for which the metallicity has been 
derived with high-resolution spectroscopy (the three clusters are shown with 
larger symbols); filled, light-blue squares indicate OCs for which [Fe/H] has 
been determined by our group. 



in their old ages an d in their large Roc (see e.g., the sample map 
in lBragagliall2007h . Be 66 is the outermost cluster in the second 
quadrant, while Be 20 and To 2 fill the large gap in Rgc between 
Be 29 (at 20.8 kpc) and Be 22 (at 14.0 kpc). Notice that in the 
literature very few cl usters have been stu died beyond a radius of 
16 kpc. For instance. ICarraro et al .1 d2007t) list only Be 31, Be 73, 
Be 25, Saurer 1, and Be 29. 

Different interpretations of the two-slopes metallicity dis- 
tribution are possible, either a normal o utcome of di s c for- 
mation and chemical e nrichment (e.g., ICarraro et al. I 120071 ; 
ISchonrich & Binnevll2009h or satellite accretion in the outer disc 



(e.g., Yong et al. 20051) . The former suggestion likely foresees a 



continuous radial distribution (even if with different slopes, maybe 
even zero, in different regions), while the latter implies some inho- 
mogeneity reflecting the satellite impact trajectory. It is then impor- 
tant to have as many clusters as possible outside the critical radius 
of 12 kpc, and to have them in different quadrants. 

We will use the information obtained here, together with the 
detailed chemical abundances from existing and future high reso- 
lution spectroscopic data, to increase the BOCCE set and derive 
conclusions based on a homogeneous analysis. 
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